Waveguide QED: Many-Body Bound State Effects in Coherent and Fock State 
Scattering from a Two-Level System 
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Strong coupling between a two-level system (TLS) and bosonic modes produces dramatic quantum optics 
effects. We consider a one-dimensional continuum of bosons coupled to a single localized TLS, a system 
which may be realized in a variety of plasmonic, photonic, or electronic contexts. We present the exact many- 
body scattering eigenstate obtained by imposing open boundary conditions. Multi-photon hound states appear 
in the scattering of two or more photons due to the coupling between the photons and the TLS. Such bound 
states are shown to have a large effect on scattering of both Fock and coherent state wavepackets, especially in 
the intermediate coupling strength regime. We compare the statistics of the transmitted light with a coherent 
state having the same mean photon number: as the interaction strength increases, the one-photon probability 
is suppressed rapidly, and the two- and three-photon probabilities are greatly enhanced due to the many-body 
bound states. This results in non-Poissonian light. 

PACS numbers: 03.65.Nk,78.67.Uh,42.50.Ct,42.50.Gy 



I. INTRODUCTION 



Recently, there has been increasing interest in designing 
quantum optical elements based on the strong coupling be- 
tween light and matter |[j]-l9|] • The strong coupling regime has 
been realized in the classic cavity quantum electrodynamics 
(QED) systems EMI, 

as well as in circuit-QED experi- 
ments 11131 — Tl6f1 . Several experimental systems have been pro- 
posed for realizing devices such as a single-photon transistor 
flU] or a quantum switch Jf], 0, [nth including surface plas- 
mons coupled to a single two-level emitter [4], a supercon- 
ducting transmission line resonator coupled to a local super- 
conducting charge qubit [6, 7], and propagating photons in a 
ID waveguide coupled to a two-level system jTa, [l9ll . Most 
of the theoretical work focuses on a single-photon coupled to 
a local quantum system modeled as a two-level system (TLS). 
The key property used in the device proposals is that, if the 
energy of the incident photon is tuned to be on resonance with 
the TLS, the system will block the transmission of photons 
due to destructive interference between the directly transmit- 
ted photon and the photon reemitted by the impurity fl, IH] . 

A more challenging task is to study the two or more photon 
scattering problem in such systems. The two-photon prob- 
lem has been addressed by Shen and Fan using a generalized 
Bethe-ansatz lfl8l [l9tl . They showed that two-photon bound 
states emerge as the photons interact with the two-level sys- 
tem. Effective attractive and repulsive interactions can be in- 
duced depending on the energy of the photons Il8ll . Such 
effective interactions between photons may provide new av- 
enues for controlling photon entanglement [20]. However, 
the scattering eigenstates were not constructed explicitly in 
Ref.[19j: the bound states were found by first constructing 
Bethe-type scattering eigenstates and then deducing the bound 
states via the completeness of the basis. It is difficult to gener- 
alize the method in Ref. 19 to solve the three-photon (or more) 
scattering problem in which we expect more complicated and 
interesting photon correlations. 
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FIG. 1: (color online) Sketch of the structure considered: a ID con- 
tinuum of bosons coupled to a two-level-system. 



Here, we present a method to explicitly construct exact n- 
photon scattering eigenstates and then use the eigenstates to 
analyze the scattering of Fock- and coherent-state wavepack- 
ets. The system consists of a ID bosonic continuum cou- 
pled to a local two-level-system as shown in FigureQ] First, 
we explicitly construct the n-photon (n — 1 to 4) scattering 
eigenstates by imposing open boundary conditions while re- 
quiring that the incoming wavefunctions consist entirely of 
plane waves I2H I22I1 . In addition to two-photon bound states, 
three-photon bound states appear in the three-photon scatter- 
ing eigenstates, and likewise n-photon bound states appear 
in the n-photon scattering eigenstates. Second, to show the 
significance of these bound states in the scattering of prac- 
tical light sources, we study the scattering of one-, two-, and 
three-photon Fock state wavepackets. It is shown that the two- 
and three-photon bound states dramatically enhance the trans- 
mission of two- and three-photon wavepackets, respectively. 
Third, we study the scattering of coherent states to determine 
the impact of the bound states on both the photon correla- 
tion and the statistics of the transmitted and reflected pho- 
tons. Strong bunching and antibunching effects appear, and 
the statistics are non-Poissonian. 

The paper is organized as follows. In Sec. II, we intro- 
duce the model, solve for the «-photon scattering eigenstates 
in the n = 1 to 4 cases, and construct the correspon ding S- 
matrix based on the Lippmann-Schwinger formalism 12311 . In 
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Sec. Ill, the impact of bound states on the photon transmission 
is studied for initial Fock-state wavepackets with photon num- 
ber of one, two, and three. In Sec. IV, we present the analysis 
of photon correlation and statistics for coherent state scatter- 
ing. Finally, we conclude in Sec. V. 



II. SCATTERING EIGENSTATES 

The system we study consists of a two-level system cou- 
pled to photons propagating in both directions in a one- 
dimensional waveguide 

HEIDI]. The system is modeled 

by the Hamiltonian lfl8ll 

/l r + d + d l 

dx—laJx) — cir{x) - al(x) — ai(x)\ 
i l K dx L dx 1 

+(e- — )|e><e| + I dxV6(x)[[a f R (x) + a f L (x)]S~ + h.c], (^) 

where a R (x) / c^ L (x) is the creation operator for a right- 
going/left-going photon at position x, e is the level splitting 
between the ground state \g) and the excited state \e) of the 
two level system, F' is the decay rate into channels other 
than the ID continuum, V is the frequncy-independent cou- 
pling strength, and S~ - \g)(e\ is the atomic lowering oper- 
ator. Throughout the paper, we set the group velocity c and 
Plank's constant H to 1 for simplicity. 

It is natural to transform to modes which are either even 
or odd about the origin, a^ a {x) = \a R {x) ± a^(-x)]/ V2. The 
Hamiltonian ([TJi is then decomposed into two decoupled 
modes: H — H„ + H„ with 



a e {x) + (e-iT'l2)\e){e\ 



H„ = 



r l t d 

I dx—a e (x)— 
J i dx 

+ J dxV S(x)[al(x)S ~ +h.c], 

r i t d 

I dx—a (x) — a (x), 
J i dx 



(2a) 
(2b) 



where the effective coupling strength becomes V = V2V. 
Note that the odd mode is free. The number operator for 
even bosons is n e - J dxa\(x)a e (x), that for the odd bosons is 

n„ - J dxa] ) (x)a n (x), and the occupation number of the two- 
level system is n t i s = \e){e\. Because H commutes with certain 
number operators, [H, n e + n t i s ] = [H, n ] = 0, the total number 
of excitations in both the even and odd spaces are separately 
conserved. We will now focus on finding the non-trivial even- 
mode solution and then transform back to the left/right repre- 
sentation. 

A n-excitation state (n-n e + n t i s ) is given by 



\'P, 



t J* 1 ...^ ie ,(, 1 ,..,„_ 1 ,4ta)-4(^ 1 )io.< ! >/3) 



where \0,g) is the zero photon state with the atom in 
the ground state. From H e \ip„) = E n \ifr„) , we obtain the 
Schrodinger equations 



[\(di + ■ ■ ■ + d„) - E n ]g n (xu - ■ ■ ,x n ) 

Vr I 

+ — \5{x\)e n (x 2 ,--- ,x„) + --- + 6(x n )e„(xi,--- ,x„-i)\ = 0, 

[t(<3i + ■ ■ ■ + d„-i) - E„ + e- iY' I2\e n (x\, - ■ ■ ,x„_i) 



+ nVg n (Q,xu--- ,x n -i) = 0, 



(4) 



where the eigenvalue E n = k\ +k 2 -i — k„, and g n (xi,--- ,x„) 
is discontinuous at jCf = 0, i = l,---n. In all the following 
calculations, we set g„(0 xu--- ,x n -\) = [g„(0 + ,x\,- ■ ■ ,x„_i) + 
g„(Q~ ,x\,--- ,x„-\)]/2 i2ll l22ll . The scattering eigenstates 
g n (xi,- ■ ■ ,x„) and e n (x\,--- ,x n -\) are constructed by im- 
posing the boundary condition that, in the incident region, 
g n {x\,--- ,x n ) is free-bosonic plane wave. That is to say, for 

gn(xi,--- ,x„) = —\h k (x Q )---h kn (x Qn ), (5a) 



h k {x) = 



V2^ 



(5b) 



For n = 1, plane-wave solutions are sufficient to satisfy 
Eq. (H]i with eigenenergy E -k: 



gi(x) = gk(x)=h k (x)[8(-x)+t k 0(x)], 

_ _ k-e+iT'/2-iT c /2 
h ~ k-e+iV l2 + iY c l2' 



(6a) 
(6b) 

(6c) 



where 8(x) is the step function and T c — V 2 = 2V 2 is the spon- 
taneous emission rate from the two-level-system to the ID 
continuum. Note that t k is the transmission coefficient for the 
even problem; because the even mode is chiral, |4| = 1 when 

r = o. 

For n — 2, plane-wave solutions are not sufficient to satisfy 
Eq. ©. As discussed by Shen and Fan lfl8i[T9ll . a two-photon 
bound state must be included to guarantee the completeness of 
the basis. Here, instead of extracting the bound state through 
a completeness check 111 <8l Il9ll . we construct the scattering 
eigenstate explicitly and find a two-photon bound state con- 
tribution to the solution, as has been done in the open inter- 
acting resonant-level model J2lll . We require the two-photon 
solution to satisfy Eq. d5al ) in the region X\,x 2 < and solve 
for the solution in other regions using Eq. (0]). This method 
of constructing scattering eigenstates can be generalized to 
three-, four-, and even more photon cases. In the Appendix, 
it is shown that the two-photon eigenstate with eigenenergy 
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E = ki + k 2 is 



<?20) 



B kp v kp 2 ( x Qi> x Q2) 



g2(x\,x 2 ) = gk u k 2 {x\,x 2 ) = —\/_ i Sk { { x Q l )8k 2 (XQ 2 ) 



= ^te(0 + ,x)-g 2 (0-^)], (7b) 

= - l)(4p 2 - 1)%, (xQ 2 )h kpi (x Ql ) 

xe { ~ r/2 - ie)]x ^~ x s l] e(xQ 2 -x Ql ). (7c) 



Here, P = (P\,P 2 ) and Q = (gi, 62) are permutations of (1,2) 
needed to account for the bosonic symmetry of the wave- 
funtion, and T -Y c + V is the total spontaneous emission 

(2) 

rate. The two-body bound-state term B k , (xq\ , xq 2 )0(xq 1 ) 



is generated when there are two photons interacting with the 
same two-level-system (TLS), while the TLS can only absorb 
one photon at one time. The binding strength of the two pho- 
tons depends on the total spontaneous emission rate T. Con- 
ceptually, two photons have two ways of going through the 
TLS. One way is to pass by the TLS independently as plane 
waves and gain a phase factor, which is described by the first 
term of g2(xi,x 2 ). The other way is to bind together and form 
a bound state, which is described by the second term. The 
formation of the bound state can be viewed as a result of stim- 
ulated emission: the first photon excites the TLS and the pass- 
ing of the second photon stimulates emission of the first pho- 
ton into the same right-going state, hence producing the bound 
state. 



For « = 3, a procedure similar to that used to solve the n ■ 
case yields 



g^{x\,X2,X3) = gk x ,k 2 M (xy,x 2 ,X3) (8a) 
= Tj gk > )#*2 ( x Qi )8k 3 (XQ 3 ) + Yj Sk Pi ( x Qi )B { k 2) ^ ^ (xq 2 , x Qi )6(xq 2 ) + J] flg ^ ^ (x Ql , x Ql , x Qi )6(x Ql )] , 



PQ 



PQ 



3/ 



e 3 (xi,x 2 ) = — — [g 3 (0 ,x u x 2 )-g 3 {0 ,Xi,x 2 )], 
V2V 



(8b) 



B 



.(3) 

k P{ ,kp 2 ,k P . 



(x Ql ,XQ 2 ,XQ 3 )=2(t kp - l)(t kp -l)(t kp - l)h kp (x Q2 )h kp (xQ 3 )h kp (xQ 3 )e<- r/2 k)]x ^ x ^6{xQ 32 )e(x Q2l ), (8c) 



r 



where P = (P{,P 2 ,Pj) and Q = (Q\,Q 2 , Q3) are permutations 
of (1,2,3) and G(x Qij ) = 6(x Qi ) - 0(x Qj ) for short. In addi- 
tion to the two-photon bound state, there emerges a three- 



body bound state B 



(3) 

kp x ,kp 2 ,k Pj 



(xQ i ,XQ 2 ,XQ 3 )9(xQ l ) in the re- 



gion x\,x 2 ,xt, > 0. Conceptually, there are three ways for 



the three photons to pass by the atom: (i) all three photons 
propagate as independent plane waves; (ii) two photons form 
a two-body bound state, while the other one propagates in- 
dependently as a plane wave; and (iii) all three photons bind 
together and form a three-body bound state. These three pro- 
cesses are described by the first, second, and third terms of 
gi{x\,x 2 , x$), respectively. 





FIG. 2: (color online) Schematic of different processes in four- 
photon scattering by a two-level system. The plane waves are repre- 
sented by wiggly lines, while the many-body bound states are repre- 
sented by the ovals. 



This simple picture can be applied to a general n-photon 
scattering process. For example, in the case of four-photon 
scattering, there are five ways for the four photons to pass 
by the atom as illustrated in Figure|2] (i) all four propagate 
as independent plane waves; (ii) two photons form a two- 
body bound state, while the other two propagate indepen- 
dently as plane waves; (iii) three photons form a three-body 
bound state, while the other one propagate independently as a 
plane wave; (iv) four photons form two independent two-body 
bound states; and (v) four photons form a four-body bound 
state. These five processes can be identified as the five terms 
of gn(x\,x 2 , Xi,X4) in the four-photon solution, which is given 
by 
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g4(xi,x2,x3,x 4 ) = ^[^^ 1 (^ei)^2( x e 2 )^3(^e 3 )^4( x e4) + X^i^ ei ^2^e2^£ ) k P Oav-^eJ^es) ( 9a ) 
q pq v 4 

+ Yj gk ^ ( x QJ B Tp„Mp,,kJ x Qi' x Q^ x Q^ x Q2) + Tj B( £, >*p, ( x Q^ x Q2^ B f P \ P (xQ 3 ,x Q4 )e(x Ql )0(xQ 3 ) 



PQ 

z 

PQ 



B 



.(4) 

»*i> 2 .^3 ."/^ 



, (x ei ,x e2 ,x e3 ,x e4 )6'(x ei )], 



e 4 (x\,X2,Xi) ■ 



4i 

V2V 



[g4(0 + ,Xi,X2,X 3 )-g4(0 ,Xi,X 2 ,X 3 )], 



(9b) 



? (4) 

k P{ ,k/> 2 ,kp 3 ,k P4 K 



BT' ,,_ ,._ ,._ (x ei ,x e2 ,xg 3 ,x e4 )s -2 (tk Pi -l)(tk P2 -l)(tk P3 -l)(tk P4 -l)hk Pl (xQ 2 )h kp2 (xQ 3 )fa 



xe( -r/2-,c)\x Q4 -x Ql \ 6{x& _ XQ3)0(xQi _ XQ2 )0( XQl - x Ql ). 

I 



(9c) 



The scattering eigenstates of a general n-photon problem 
can be constructed recursively in a similar way: the only un- 
known term in g n (x\,--- ,x n ) is the n-photon bound state as 
all the other terms can be constructed from the solutions of 
the 1,2,- - ,(n- l)-photon problems. We extrapolate from the 
results of n — 2-4 that, for general n (> 2), the n-body bound 
state assumes the form 



n-l 



B k . 



xh kl (x n )h k2 (x 2 )- ■ -hk^xn-^hk^e^ 12 -^"-^. (10) 

We have verified this expression for n — 5. Thus we have given 
explicit formulas for constructing the exact n-photon scatter- 
ing eigenstates. 

The exact scattering eigenstates can be used to construct the 
scattering matrix. According to the Lippmann-Schwinger for- 
malism [23], one can read off the "in" state (before scattering) 
and the "out" state (after scattering) of a general n-photon S- 
matrix from g„(x\ , ■ ■ ■ , x n ) in the input region (x\ < 0, • • • ,x n < 
0) and in the output region (x\ > 0, ■ • ■ ,x n > 0), respectively. 
The "in" and "out" states of one and two photon scattering 
matrices are given by 



- J dxhk(x)al(x)\Q) 

- J ' dxtkhk(x)al(x)\0), 



(11a) 
(lib) 



and 



i*2>>. 



f dxidx 2 — [Y\hk l (xQ l )hk 2 (xQ^)\al(xi)al(x2)\0) 
J Q 

= J ' dxidx2^\^Jk { tk 1 hk i (x Ql )hk 2 (xQ 2 ) 



Q 

+ 2^ Bk P[ Jcp 3 (XQj , XQ 2 )]al(x l )al(x 2 )\0}, (12) 
PQ 



and similarly for three and four photons. The corresponding 
S -matrices are 



dh---dk n ^\$l t )ee(<!$\ 



(13) 



Notice that the unitarity of the S-matrix is automatically satis- 
fied since the incoming state \<p\ n ) e is a complete basis set in 
the even space Il9l 12311 . 

The S-matrix in the odd space is just the identity operator 
because the odd mode is free and decoupled from the impurity 
and the even mode, 



1 



S',; n = 1 dk r --dk n ^\<p { £) 00 {<p™\, 



(14a) 



/l " 
dxi ■ ■ ■ dx n — ^ Y [ h h (xQi)at(xi)\0) . (1 4b) 



Q i=l 



Finally, we wish to construct the scattering matrix in 
the right/left representation based on the S-matrices in the 
even/odd representation. For a general n-photon scattering 
problem, the possible scattering channels are that i photons 
undergo scattering in the even space and n - i photons undergo 
scattering in the odd space, with i running from to n. In ad- 
dition, the even and odd spaces are decoupled from each other. 
Therefore, the n-photon S-matrix is 



(15) 



We will use this S-matrix to study the scattering of Fock states 
and coherent state wave packets in the right/left space in the 
subsequent sections. 



III. SCATTERING OF FOCK STATES 

In order to show the significance of the many-body bound 
states, we study the scattering of a Fock state off of a two- 
level system. We assume that the incident mode propagates 
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to the right and the two level system is initially in the ground 
state. We use the S-matrices defined in Eq. ( TT3T > to evaluate 
the transmission and reflection coefficients. In practice, any 
state that contains a finite number of photons must have the 
form of a wave-packet. Thus, we start with the definition of 
the continuous-mode photon wave-packet creation operator in 
momentum space B24I1 

at = J dka(k)a\k)\0), (16) 

with the normalization condition fdk\a(k)\ 2 = 1. The corre- 
sponding continuous-mode n-photon Fock state is 




0.2 0.4 0.6 0.8 1 



\n a ) = ^|0> , 
and the output state after it scatters off the TLS is 



ryj\ FIG. 3: (color online) Single-photon transmission (P R ) and reflec- 
tion (P L ) probabilities as a function of coupling strength V. The 
incident photon is on resonance with the two level system (ko = e) 
and we have considered the lossless case T' = 0. (A = 0.1.) 



|out^>=S (n) K>. 



(18) 



A. Single-Photon Fock State Scattering 



To obtain the scattering probabilities of a Fock state 
from the S-matrix found in Section II, we follow the fol- 
lowing general procedure, (i) First, we write an n-photon 
input Fock state traveling to the right in momentum 
space: \n a ) - (1 / yfn\) J dki ■ • • dk„ a(k\)---a{k n )\k\,--- ,k„); 
(ii) Next, we apply the S-matrix on the input 
state and find the output state |out^'') = S M \n a ) - 
(1/Vn!) fdk l ---dk n a(k l )---a(k„)S M \k u --- ,k n ) in 
the even/odd basis; (iii) We transform back to the 
right/left basis. Then we project the output state 
onto the n-photon (right/left-going) momentum basis 
\h,--- ,k„) R ,--- ,\h,--- ,ki) R ® \ki+u--- ,k n } L ,--- ,\h,--- ,k„) L 
and take the absolute value square to obtain the probabilities 
P(k\,---,k n ) of finding the output state in \ki,- ■ ■ ,k n ); (iv) 
Finally, we integrate P(k\,--- ,k„) over k\,--- ,k n to obtain 
the total transmission and reflection probabilities. Here, 
a right/left-going state is defined by a positive/negative 
momentum, i.e., k\ > 0, ■■■ ,k n > for \k\,--- ,k n )R and 
ky <(),••• ,k„ <0for \k\,--- ,k„) L . 

For convenience, we choose Gaussian type wavepackets 
with the spectral amplitude 



The probabilities of transmission (P^ ) and reflection (P^P) 



for a single-photon Fock state are found as 

= f dk \(k\out ( a } )\ 2 = f dkm 

Jk>0 Jk>0 



(kf\t k \ 2 , (20a) 



^ = f dk\(k\ont^}\ 2 ^ f 

Jk<0 Jlc 



dka(kf\r k \ 2 , (20b) 



k>() 



where t k = (h + l)/2 and r k = (tk - l)/2 and f]t is the transmis- 
sion coefficient defined above for the even mode [Eq. d6cbl. 

Note that the propagation of a single-photon is strongly 
modulated by the TLS as we turn on the coupling. In the 
strong-coupling limit, a single-photon is perfectly reflected 
and the two-level atom acts as a mirror. This perfect reflection 
is due to destructive interference between the directly trans- 
mitted state and the state re-emitted from the TLS. A single- 
photon transistor [4] and a quantum switch J6J] have been pro- 
posed based on this perfect reflection. 



B. Two-Photon Fock State Scattering 



a W = (2,rA 2 )- 1/4 exp(-^^). (i 9 ) 

For all of the numerical examples in this paper, we choose 
ko = e: the central frequency of the wavepacket is on reso- 
nance with the TLS, a condition which makes the interaction 
between the photons and the TLS strongest. We take the cen- 
tral momentum ko » A so that the narrow-band condition is 
satisfied. In particular, we choose A = 0. 1 . However, we em- 
phasize that all the conclusions we draw are independent of 
the choice of A. That is because all the transmission and re- 
flection probabilities are functions of F/A, where F = 2V 2 . A 
different choice of A does not change any of the qualitative 
results, but merely rescales the spontaneous emission rate F. 



For two incident photons, following the general procedure 
above, we find that the transmission and reflection probabili- 
ties are 



3 (2) 
RR 



f ■ 

Jki>0,k 2 >0 



1 



dkidk 2 —\(kuk 2 \out {2) )\ 2 , (21a) 



RL 



-J 

J/q 

-J 

Jk v 



>0,it 2 <0 



2! 



dkidk 2 \(ki,k 2 \out (2) )\ 2 , (21b) 



3 g = dhdk 2 -\{k u k 2 \0Mi (2) )\ 2 , (21c) 



(2) f2) (2.) 

where P^, and P^[ are, respectively, the probability for 
two photons to be transmitted (right-going), one transmitted 
and one reflected, and two photons reflected (left-going). 
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To show the significance of the bound state in the propaga- 
tion of multi-photon Fock states, we separate each of the prob- 

(2) (2) (2) 

abilities P RR , P\l and P\l into two parts. One part is the con- 
tribution from only the plane wave term (labeled PW), which 
is the direct transmission or reflection. The other is the con- 
tribution from all the other terms (labeled BS), including the 
bound state term as well as the interference term between the 
plane wave and bound state. Notice that the BS part vanishes 
in the absence of bound state, as in the case of single-photon 
scattering. Therefore, it is a manifestation of the nonlinear 
effect caused by the interaction between the TLS and two or 

(2) 

more photons. As an example, P ~ 
is 



RR 



split into PW and BS parts 



D<2) 

RR 



(Opw 



t(h,k 2 ) 
B(k u k 2 ) 



-I ■ 

Jki>0,k 2 >0 



dk l dk 2 \t(k l ,k2) + B(kuk 2 )\ 2 (22a) 



Opw + (Obs> 



-L • 

Jki>0,k 2 >0 



dk\dk 2 \t(k\,k 2 )\ , 



(22b) 
(22c) 



( dkidk 2 [t*(k u k 2 )B(k u k 2 ) (22d) 



- I 



+t(kuk 2 )B*(k u k 2 ) + \B(k u k 2 )\ 2 
a(ki)a(k 2 )t kl t k2 , 

—i/2jr —i/2n n 



+ ■ 



(22e) 
(22f) 



ki-e+'-j k 2 -e+'-j 

x I dk a(k )a(ki +k 2 -k )r k > r, . _ k > . 
Jk'>o 

Figure|4] shows the three transmission probabilities p£L 

(2) (2) 

P R jj and P\l for our standard parameters, with the contribu- 
tions from the plane wave and bound state plotted separately 
in panels (a)-(c). Note that the presence of the bound state has 
a very substantial effect on these transmission probabilities. 

(2) (2) 

As shown in panels (a) and (b), P RR and P R ^ are enhanced 
by the formation of the bound state. This is mainly due to 
constructive interference between the plane wave and bound 

(2) 

state. In contrast, panel (c) shows that P L ^ is strongly reduced 
in the presence of the bound state because of destructive in- 
terference between the plane wave and bound state (change 
from ~ 0.8 to ~ 0.4 at V — 0.5). Therefore, the presence of the 
bound state tends to increase the one-photon and two-photon 
transmission, while suppressing the two-photon reflection. 

A particularly interesting aspect of the results in Fig|4] is 
that the effect of the bound state is most prominent in the 
intermediate coupling regime, not at the strongest coupling. 
This is because, first, in the weak coupling limit, the inter- 
action is too weak to produce a pronounced bound state for 
two-photon scattering, while, second, in the strong coupling 
limit, the TLS responds to the first photon too quickly (in a 
duration of order 1 /F with T = 2V 2 ) for the second photon 
to produce a significant nonlinear effect. (The formation of 
the bound state requires the presence of both photons at the 
two-level system.) The optimal coupling strength V„, for pro- 
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FIG. 4: (color online) Two-photon transmission and reflection proba- 
bilities as a function of coupling strength V. (a) Probability that both 

(2) 

photons are transmitted (and hence are right-going, P^). (b) Prob- 
ability that one photon is transmitted and one reflected (right-left, 

(2) 

P RL )- (c) Probability that both photons are reflected (both left-going, 

(2) 

lD' W The three processes on a single plot. The label PW refers 
to the contribution from the plane- wave term only, while BS refers 
to all the other contributions involving bound-state terms. The inci- 
dent photons are on resonance with the two-level system (kg = e), we 
consider the lossless case T' = 0, and A = 0.1. Notice the large effect 
of the bound state on these quantities. 



ducing nonlinear (bound state) effects lies at intermediate cou- 
pling, when the spontaneous emission rate T is on the order of 
the wavepacket width A CV m ~ 0.4 when A = 0.1). 



C. Three-Photon Fock State Scattering 

Following the general procedure for obtaining scattering 
probabilities, the transmission and reflection probabilities for 
three-photon Fock state scattering are defined as 



,(3) 
RRR 

(3) 
RRL 

(3) 
RLL 

(3) 
LLL 



j 

Jk[> 

J 

r 

•Jki> 

J 



1 



dkidk 2 dkj, ^Kh^Mloat^)] 2 , 



0,k 2 >Q,k 3 >0 



3! 
1 



(3)\i2 



0,/t 2 >(U-3<0 



dk\dk 2 dk^ —\(^\,k 2 ,k^\ovX a )\ 



2! 
1 



dk\dk 2 dh ^7l<^i,^2,^3|out^ 3) >| 2 , 



0,/t 2 <0,A- 3 <0 



2! 
1 



dk\dk 2 dk 3 --|<Ari, A:2,^3|out^ 3, )| 2 , 



<0,it 2 <0,<: 3 <0 



3! 



(23) 



(3) (3) (3) (3) 

where P RRR , Prrl' rlv anc ' lll me tne probabilities for 
three photons being transmitted (all right-going), two trans- 
mitted and one reflected, one transmitted and two reflected, 
and all three reflected (left-going), respectively. As in the two- 
photon scattering case, we separate each probability into two 
parts: the contribution of only the plane wave term (labeled 
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FIG. 5: (color online) Three-photon transmission and reflection probabilities as a function of coupling strength V. (a) Probability of all three 

(3) (3) 
photons transmitted (P RRR )- (b) Probability of two photons transmitted and one reflected (P RRL )- (c) Probability of one photon transmitted and 

two photons reflected (P R l L )- (d) Probability of all three photons reflected (P^l L )- (e) all together. The label PW refers to the contribution 
from only the plane wave term, while BS refers to all the other contributions, involving bound state terms. The incident photons are on 
resonance with the two level system (jfco = e), we consider the lossless case T' = 0, and A = 0. 1. Note the large bound state effects. 



PW) and the contribution from all the other terms (labeled 
BS), including the bound states as well as the interference be- 
tween the plane wave and bound states. The probabilities and 
the decomposition into PW and BS parts are plotted in Fig- 
ure|5]for our usual parameters. 

Figure|5] shows that the bound state contribution to the 
transmission probabilities is, as for two photons, very substan- 
tial. In panels (a) and (b), the BS parts of P RRR and P RRL are 
positive; thus, these probabilities are enhanced by the bound 

(3) 

states. Panel (d) shows that P L } L is suppressed by the bound 
state contribution for arbitrary coupling strength. In contrast, 

f31 

as we increase the coupling strength, P r [t is first suppressed 
and then enhanced by the BS part as shown in Figure|3c). 
Tuning the coupling strength changes the relative phase be- 
tween the plane wave and bound state parts; for P R [ L , the in- 
terference between them happens to change from destructive 
to constructive as the coupling strength increases. Finally, as 
in the two-photon case, the most pronounced bound state ef- 
fects occur in the intermediate coupling regime instead of the 
strong coupling limit. 

To sum up this section, we point out that all the curves 
plotted in Figs. [3]|5] are universal in terms of the choice of A. 
Because A appears in the scattering probabilities [P R \ etc.] 
only in the ratio T/A, a different choice of A (i.e., other than 
0.1 used in the figures) is equivalent to rescaling V and does 
not change the shape of the curves. Therefore, the substan- 
tial bound state effects observed here are intrinsic for multi- 
photon scattering processes in this system, independent of the 
details of the wavepackets. 



IV. SCATTERING OF COHERENT-STATES 

We now turn to studying the scattering of coherent states 
in order to show, first, the strong photon-photon correla- 
tion induced by the the two-level system and, second, the 
change in photon number statistics. The incident coherent 
state wavepacket is defined by II24I1 



with a\ = j dka(k)a'(k)\0), and mean photon number n 
J dk\a(k)\ 2 . A Gaussian type wavepacket is chosen 



a(k)- 



V» / (k~ko) 2 \ 
(2^A 2 )!/ 4eXP ' 4A 2 ' 



(25) 



for numerical evaluations, we use, as before, A = 0.1 and ko ■ 
e » A. The output state lout,*} is then 



iout ff > = £V"V> 



(26) 



We assume the incident coherent state is right-going and the 
two-level system is in the ground state initially. We present 
the analysis of second-order correlation and photon number 
statistics in the transmitted field. 



A. Correlation 

The second-order correlation function of the transmitted 
field is defined as [24 



(2) (out a \al(x l )al(x 2 )a R (x2)a R (xi)\outa) 

g i R \x2-x l ) = 2 1 . (27) 

<out ff |aj,(.xi)a/;(.ri)|out a > 2 

We consider the mean photon number h < 1.0. In this case, 
the probability to find n > 3 number states is much smaller 
than that of n — 2 number states. Moreover, the contributions 
from n > 3 states to are at least one order of A (= 0.1) 
smaller than that from the n — 2 state. Therefore, we neglect 
the contributions from n > 3 number states. The second-order 
correlation function simplifies to 



\j dkidk2a(ki)a(k 2 )(tk 1 tk2 



\J dk\dkia{k\)a{ki)t kx t kl 



\a) = e a --" /2 \0} 



(24) 



(28) 

The contributions from the directly transmitted state and 
the bound state can be identified as the first term and sec- 
ond term in the numerator of g R \x2 - x\) in Eq. d28l>. In 

(2) 

the absence of the bound state, g R (x 2 — JCj) is always equal 
to unity. As we turn on the interaction, the interference be- 
tween the directly transmitted state and the bound state will 
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FIG. 6: (color online) Second-order correlation of the transmitted 
field given an incident coherent state with n < 1 at various cou- 
pling strengths V to the ID continuum, (a) V = 0.16, (b) V = 0.26, 
(c) V = 0.34, (d) V = 0.38, (e) V = 0.40, (f) V = 0.45. The sponta- 
neous emission rate to channels other than the ID continuum is set 
to T' = 0.10. Notice that the correlation behavior is very sensitive to 
the coupling strength to ID continuum, showing both bunching and 
antibunching. 



give rise to interesting correlation behavior. Figure|6] shows 
the second-order correlation as a function of T{x2 - x\) at var- 
ious coupling strengths, V, to the ID mode with P = 0.1. In 
the weak coupling limit (V = 0.16) as shown in Figure|6ja), 

(2) 

the directly transmitted state dominates and g~ (0) is slightly 
smaller than 1. We observe a slight initial antibunching. As V 
increases [Figure|6jb)-(c)], g^(0) further decreases and the 
initial antibunching gets stronger and becomes strongest at 

(2) 

V — 0.34 when g R (0) — 0. Notice that the antibunching is get- 
ting weaker as one moves away from the origin for V < 0.34. 
Further increase of V starts to change the initial antibunching 
[V = 0.38, gf(0) < 1] to bunching [V = 0.45, gf(0) > 1] as 
shown in Figure[6|d)-(f). In this case, the bound state starts 
to dominate the correlation behavior. It is remarkable that, 
for V > 0.34, the initial antibunching (V < 0.40) or bunch- 

(2) 

ing (V > 0.40) is followed by a later antibunching g R (0) = 0, 
which is caused by the cancellation of the directly transmitted 
state and the bound state. The formation of the bound state 
gives rise to a rich phenomenon of photon-photon correlation, 
which is very sensitive to the coupling strength V to the ID 
mode. Effective attractive or repulsive interaction between 
photons is induced by the presence of a single two-level sys- 
tem fill . 

Our findings agree with the results obtained by Chang et al. 
10] using a very different approach. In the lossless F' = case, 
as we increase the coupling strength, the transmission for indi- 
vidual photons is reduced rapidly [see, for example, Figure[3] 
and Figure|4ja)]. But the two-photon bound state can strongly 
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FIG. 7: (color online) Photon number distribution of the transmitted 
field compared with a coherent state. We considered the lossless 
case T' = 0. The statistics is non-Possonian with the 2 and 3 photon 
content enhanced. 



enhance the transmission. Therefore, we will observe a strong 
initial bunching followed by a later antibunching, similar to 
FigurelSff). 



B. Photon Number Distribution 

Given the output state |out ff ), we measure the photon num- 
ber distribution in the transmitted field following the general 
procedure described in Sec. III. 



P{) 
Pi 

Pi 

Pi 



= |<out ff |(|0) s ®|/) L )| 2 , 



f 
f 

Jk\J(2> 

f 

Jki,k 2 ,k 



dk\(0Ut a \(\k)R®\I)L)\\ 



1 9 

dhdki—KoutalQkiMn®]!)!)], 



'2! 



1 



dk l dk 2 dk 3 — \{out a \(\kuk2,k 3 ) R ®\I) L )\ z , 



3! 



(29) 



where is the complete basis set in the left-going photon 
space. We consider a mean photon number h < 1.0 in the in- 
cident coherent state. In this case, the probability to find the 
four photon state is negligible (< 1 .6%). We compare the pho- 
ton number distribution P n of the output state with (P ; ,)p isson 
of a coherent state having the same mean photon number. 

Figure|7] shows the ratio between (P„)p isson and P n as a 
function of the coupling strength V and the mean photon num- 
ber n of the incident coherent state. The zero-photon proba- 
bility does not deviate from that of a coherent state much in 
the whole parameter region we considered. The one-photon 
probability is smaller than the corresponding probability in a 
coherent state. In contrast, the two- and three-photon proba- 
bilities are much larger than the ones in a coherent state, espe- 
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cially in the strong coupling regime. This is to say, the inter- 
action between photons and the two-level system redistributes 
the probabilities among different photon numbers. The one- 
photon probability is reduced and is redistributed to the two- 
and three-photon probabilities. This is mainly because the 
bound states enhance the transmission of multi-photon states 
as we have shown in Sec. Ill B and C. In conclusion, we ob- 
tain a non-Poissonian light source after the scattering. It is 
perhaps possible to use this strongly-correlated light source to 
perform passive decoy-state quantum key distribution in order 
to raise the key generation rate 025142811 . 



V. CONCLUSION 



In this paper, we present a general method to construct the 
exact scattering eigenstates for the problem of n-photons in- 
teracting with a two-level system. Many-body bound states 
appear in the presence of the coupling between photons and 
the two-level system. Furthermore, the scattering matrices are 
extracted using the Lippmann-Schwinger formalism. We em- 
phasize that the completeness of the S-matrices is guaranteed 
by imposing open boundary conditions and requiring the inci- 
dent field to be free plane waves. Based on the S-matrices, we 
study the scattering of the Fock states and coherent states. The 
bound states are shown to enhance the transmission of multi- 
photon states and suppress the transmission of single-photon 
states. In the transmitted field of coherent state scattering, the 
photons exhibit strong bunching or antibunching effects de- 
pending on the coupling strength. This is a manifestation of 
the many-body bound states. Finally, we determine the pho- 
ton number distribution and find that the one-photon state is 
transferred to two- and three-photon states. This results in a 
non-Poissonian light source which might have applications in 
quantum information. 
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Appendix: Two-photon Scattering Eigenstate 

In this appendix, we show in detail how we obtain the two- 
photon scattering eigenstate [Eq. (0] by imposing the open 
boundary condition Eq. (0. The equations of motion for the 



two-photon case read 

[t(3i + d 2 )-E 2 ]g 2 (xi,x 2 ) + 



V 



-[6(xi)e2(x2) + 6(x2)e 2 (xi)]=0, (Ala) 



[- — -E 2 + e- ir/2]e 2 (x) + 2Vg 2 (0, x) = 0, (Alb) 



i dx 



which can be cast into the following set of equations 

[j(d 1 + d 2 )-E 2 ]g 2 (x 1 ,X2) = 0, (A2a) 

(A2b) 



2i 



e 2 (x) = -^[g 2 (0 + ,x)-g 2 (0-,x)], 

[~-E 2 + e-ir/2]e 2 (x) 

+V[g 2 (0 + ,x)+g 2 (0-,x)]=0, (A2c) 
e 2 (0 + ) = e 2 (0-). (A2d) 

Here, g 2 (x\, X2) is discontinuous at x\ — Q,X2 = and we set 
g2(x,0) = [g2(x,0 + ) + g2(x,0~)]/2. We eliminate e 2 (x) from 
the above equations and obtain 



[-(d 1 +d 2 )-E 2 ]g 2 (x 1 ,x 2 ) = 0, 



(A3a) 



[I L - e 2 + e - iT 12 - iT c l2]g 2 (Q + , x) 



i dx 



[)^- E 2 + e-iT'l2 + iT c /2]g 2 (0-,x), (A3b) 



£ 2 (o + ,o + ) - * 2 «r , o + ) = g2 (o-, o + ) - * 2 «r ,o-). (A3c) 

Because of the bosonic symmetry, we can solve for g 2 (xi , x 2 ) 
by first considering the half space x\ < X2 and then extend- 
ing the result to the full sapce. In this case, there are three 
quadrants in real space: Q x\ < X2 < 0; © x\ < < x 2 ; ® 
< x\ < x 2 , Eq. ( IA3bb can be rewritten as two separate equa- 
tions 



[~ - E 2 + e-T'/2 - iT c /2]gf(x,0 + ) 



\-—-E 2 + e-iT'/2 + ir c /2]g®(x,0-), forx<0, (A4a) 
1 i dx 1 z 



[~ ~E 2 + e- iV 12 - iT c /2]gf(0 + ,x) = 

\- — -E 2 + e-iT'/2+ir c /2]gf(0-,x), forx>0. (A4b) 
1 i dx 1 L 

Substituting g®(x\,0~) [Eq. (0] into Eq. dA4a| ), we solve to 
find 

gf(x,0+) = ±[t k2 e -^ + t ki e -^^ (A5) 
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where A is a constant to be determined. Applying the con- 
straint Eq. dA3at to g®(x,0 + ), we obtain 

I e '(k l x l +k 2 x 2 ) _ e i(k 2 x i +k l x 2 ) 
82 (XUX2) = - [k 2 ^ + t kl ^ \ 

+ A e (T/2+i£)(x 2 -x l ) e i{k\ +k 2 ) xi 

From Eq. ( IA6I 1. we can identify A to be zero: otherwise, the 
solution is not normalizable [e T ^~ x ^' 2 is divergent when X2 - 
x\ — > 00], Hence, §2(^1,^2) in region © is given by 

1 e i(kixi+k 2 X2) _ e Khx\+kxX2) 

gf(xi,x 2 ) = -[ik 2 +t kl j. (A7) 

Substituting Eq. (IA71 > into Eq. dA4bb yields 

g®(0+ , x) = If*, t k2 [^ + e -^] + Set-W> (A8) 



where B is a constant to be determined. Again, applying the 
constraint Eq. ( IA3ab to g®(0 + ,.x), we obtain 

n 1 - - [ e i< - k ' x,+k2X2 '> e Kk\x 2 +k 2 x\) 

+ Be^~ r/2 ~ i€ ^ X2 ~ X[ ) e i{kl +k ^ X2 . ( A9) 

Finally, B is found by substituting Eq. (O, Eq. ( |A7t , and 
Eq. dA9l ) into the continuity condition Eq. flA3c| >. yielding 

(ft - 1)(4, - 1) 
B = -- — -. (A10) 

Extending these solutions from the half space to the full space 
using the bosonic symmetry gives rise to the two-photon scat- 
tering eigenstate given in Eq. (0 of the main text. 
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